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0\ . ABSTRACT 

qj ', Proton-proton correlations and correlations of protons, deuterons and tritons with 

a particles from spectator decays following ^^^Au + ^^^Au collisions at 1000 MeV per nu- 
^ ' cleon have been measured with two highly efficient detector hodoscopes. The constructed 

C ' correlation functions, interpreted within the approximation of a simultaneous volume de- 

cay, indicate a moderate expansion and low breakup densities, similar to assumptions made 
^ . in statistical multifragmentation models. 



PACS numbers: 25.70.Pq, 21.65.+f, 25.70.Mn, 25.75.Gz 



Expansion is a rather basic conceptual feature of the multifragmentation of heavy nuclei. 
A volume of about three to eight times that occupied at saturation density is assumed in 
the statistical models aiming at a phase space description of the multi-fragment breakups 
im, 1^. Expansion also provides the link to the nuclear liquid-gas phase transition; only 
one third of the saturation value is expected for the critical density of nuclear matter |Q. 
The experimental confirmation of expansion to low breakup densities is therefore of the 
highest significance for the understanding and interpretation of the multifragmentation 
phenomenon [^. 

In central collisions of heavy nuclei, rapid expansion is evident from the observation of 
radial collective flow [Q . For the fragment decay of spectators following collisions at rela- 
tivistic energies, the case studied in this work, significant radial flow has not been observed 
IP]. Here evidence for expansion has been obtained, indirectly, from model comparisons. 
Models that assume sequential emission from nuclear systems at saturation density under- 
predict the fragment multiplicities while those assuming expanded breakup volumes yield 
satisfactory descriptions of the populated partition space 0, H, |, |^. The disappearance of 
the Coulomb peaks in the kinetic-energy spectra of emitted light particles and fragments, 
associated with increasing fragment production, provides additional evidence consistent 
with volume emission from expanded systems [|1^, Q. There are also other dynamical 
and statistical observables that have been interpreted as evidence for expansion in recent 
papers [0,0,111. 

Interferometry-type methods permit experimental determinations of the breakup vol- 
ume or, more precisely, of the space-time locations of the last collisions of the emitted 
products [0. In the nuclear regime, correlation functions for light charged particles, pre- 



dominantly proton-proton correlations, have been widely explored for that purpose |jT^ |TI 
Depending on the assumed reaction scenario and energy regime, both time scales and 
breakup radii have been deduced. The time scales for the decay of highly excited spectator 
nuclei produced at relativistic bombarding energies should be rather short |19[ , and we may 



expect that the correlation functions are mainly sensitive to the spatial extension of the 
source. More importantly, if we assume a rapid volume breakup of the system, the quantity 
of interest will be the local density, i.e. the mutual proximity of the nascent fragments and 
light particles. These densities are obtained in the limit of assuming a zero-lifetime in the 
source analysis. 

In this Letter, we present the results of correlation measurements for spectator decays 
following collisions of ^^^Au -|- ^^^Au at a bombarding energy of 1000 MeV per nucleon. 
Besides proton-proton coincidences, also coincidences of protons, deuterons, and tritons 
with a particles were measured and correlation functions were constructed. For their 
quantitative interpretation, it is assumed that they are dominated by the effect of final- 
state interactions. The results are found to be consistent with low breakup densities with 
values close to those assumed in the statistical multifragmentation models. 

Beams of ^^^Au with incident energy 1000 MeV per nucleon were provided by the heavy- 



ion synchrotron SIS and directed onto targets of 25-mg/cni^ areal thickness. Two muhi- 
detector hodoscopes, consisting of a total of 160 Si-CsI(Tl) telescopes in closely-packed 
geometry, were placed on opposite sides with respect to the beam axis. The angular range 
6iab from 122° to 156° was chosen with the aim of selectively detecting the products of the 
target-spectator decay. Each telescope consisted of a 300-/im Si detector with 30 x 30 mm^ 
(96 detectors) or 25 x 25 mm^ (64 detectors) active area, followed by a 6-cm long CsI(Tl) 
scintillator with photodiode readout. The distance to the target was about 60 cm. 

The products of the projectile decay were measured with the time-of-flight wall of the 
ALADIN spectrometer and the quantity Zhound was determined event-by-event. Z^yound 
is defined as the sum of the atomic numbers Zi of all projectile fragments with Zi > 2. It 
reflects the variation of the size of the primary spectator nuclei and is inversely correlated 
with its excitation energy. Because of the symmetry of the collision system, the mean 
values of Zhound ^OT the target and the projectile spectators within the same event class 
have been assumed to be identical. 

Examples of correlation functions constructed for the four types of coincidences p-p, 
p-^He, d-'^He, and t-'^He are shown in Fig. 1, sorted into four bins of Zbound as indicated. 
The uncorrelated yields were obtained with the technique of event-mixing and normalized 
in the range of large relative momenta q > 70 MeV/c (p-p) and q > 150 MeV/c (p-a, d-a, 
and t-a). In the off-line analysis, thresholds were set at Eiab > 20 MeV for all particles 
p, d, t, and a and all hodoscope detectors. The shapes of the correlation functions are 
sensitive to the threshold. For p-p, e.g., the suppression at small q tends to disappear 
if the threshold is set at Ei^b > 10 MeV or lower, presumably as a result of increasing 
contributions from evaporation and sequential decay. These long-lifetime components are 
suppressed if higher thresholds are chosen. The pairs of particles were also requested to be 
detected in the same hodoscope in order to avoid correlation effects that appear at large q 
and are beheved to be due to a collective (sidewards) motion of the proton-emitting source. 

The p-p correlation functions are characterized by a depression at small relative mo- 
mentum and by a weakly pronounced peak near relative momentum g = 20 MeV/c, caused 
by the S-wave nuclear interaction and used for the quantitative interpretation. The three 
correlation functions of the hydrogen isotopes with a particles are dominated by the res- 
onances corresponding to the ground state of ^Li and by the 2.19-MeV and 4.63-MeV 
excited states of ^Li and ^Li, respectively. The observed widths of the ^Li and ^Li peaks 
in the d-a and t-a correlation functions represent the experimental resolution which is 
mainly determined by the angular resolution following from the geometry of the detector 
hodoscopes. A striking feature of the data is the overall stability of the peak heights as a 
function of Zhound which indicates source extensions that do not change dramatically with 
impact parameter. The largest deviations of the peak height and of the overall shape of 
the correlation functions appear in the bin 60 < Zhound < 79, corresponding to the largest 
impact parameters. 

The analysis of the p-p correlation functions was performed with the Koonin-Pratt 
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Figure 1: Correlation functions constructed for p-p, p-^He, d-^He, and t-^He coincidences 
(from top to bottom) for spectator decays following ^^^Au + ^^^Au collisions at 1000 MeV 
per nucleon. The data are sorted according to Zbound into four bins with limits given in 
the top panel of each column. The lines represent the results of the calculations used to 
extract source radii (see text). Note that the scales of the abscissa are different for the top 
row (p-p correlations) and for the remaining three rows of panels. 



formalism [2D| in the zero-lifetime limit. In this form, the analysis includes the effects of 
quantum statistics and of the mutual nuclear and Coulomb final-state interactions but it 
ignores the long-range Coulomb repulsion of the two protons from the emitting source. In 
order to assess the magnitude of the latter effect, classical Coulomb trajectory calculations 
were performed. In addition, calculations with the three-body quantum model of Lednicky 
et al. ||2ll] were used to identify possible systematic uncertainties. A uniform sphere with 
radius R and statistical momentum distributions corresponding to the measured kinetic- 
energy spectra were assumed for the proton source; the model results are slightly dependent 
on the particle momenta through the applied experimental filter. 

The quality of the obtained results and the sensitivity to the radius parameter R is illus- 
trated in Fig. 2 for the case 20 < Zhound < 40. With the filtered Koonin-Pratt calculations, 
the most satisfactory description of the data is obtained with R = 8.2 fm (full line). The 
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Figure 2: Proton-proton correlation func- 
tion for the interval 20 < Zhound < 40 
in comparison with the predictions of the 
Koonin-Pratt formalism before (dotted) 
and after applying the experimental fil- 
ter (full line), and after correction for 
the Coulomb repulsion from the emitting 
source (dashed line, closely following the 
full line). The calculations were performed 
for a homogeneous spherical source with 
radius R = 8.2 fm. 



height and the width of the resonance peak are rather well reproduced. For comparison, 
the results of the same calculation before filtering (dotted line) and after correction for 
the Coulomb repulsion from the source (dashed line) are also shown. The filtering effect 
is equivalent to a change of the radius by AR ^ 0.3 fm whereas the modification caused 
by the charge of the emitting source is almost negligible. The latter is not unexpected in 
the present case, especially for the more central collisions which produce spectators of very 
moderate total charge (e.g. ~ 45 for 20 < Zhound < 40). In addition, with the assumption 
of a simultaneous volume break-up, only the charge inside the volume corresponding to the 
radial position of the particle was assumed to contribute to its acceleration. Correlation 
peaks of nearly the same height but with slightly smaller width were obtained with the 



three-body quantum model [^. In the two-body approximation, the difference between 
the results obtained with the Koonin-Pratt and the Lednicky et al. formalisms is prac- 
tically negligible. Therefore, and because of the reduced importance of the three-body 
Coulomb effect, the quantitative analysis of the measured p-p correlation functions was 
performed with the filtered Koonin-Pratt simulations. 

The best fits, generated by minimizing y^ within the range 10 < g < 35 MeV/c are 
shown in Fig. 1. The corresponding radii, listed in Table 1, are close to about 8 fm, up 
to nearly 9 fm for the bin of largest Zbound which, however, seems to be afflicted with the 
largest uncertainties. These values are distinctly but not excessively larger than the radius 
R = 6.7 fm of a gold nucleus at normal density 0.16 fm~'^. The quoted errors are purely 
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Table 1: Mean mass numbers < A > and extracted source radii in units of fm for the 
indicated bins of Zbound- The regions of relative momentum q selected for the x^ test 
are 10 - 35, 30 - 70, and 30 - 54 MeV/c for the p-p, p-a, and d-a correlation functions, 
respectively. The numbers in brackets denote the minimum x^ per degree of freedom. 



statistical and were determined from the radii for which x^ exceeds its minimum by an 
amount equal to the minimum value of x^ per degree of freedom (i.e. for x^ = Xmm'''^/('^~1) 
where n is the number of data points included). The additional systematic uncertainty, 
mainly resulting from the arbitrariness of choosing the normalization interval and estimated 
to be about 0.5 fm, is much larger. 

Correlation functions for proton pairs belonging to four intervals of the laboratory pair 
momentum Psum = |pi +P2I, integrated over Zhound, are shown in Fig. 3. Only the data 
from the 96-element hodoscope were used. The combined effects of the energy threshold 
Eiab > 20 MeV, of the binning in pair momentum, and of the finite solid-angle acceptance 
of the detector hodoscope cause the limitation of the populated q range for the bins of 
smaller Pgum- A normalization interval 70 < g < 90 MeV/c was therefore chosen for these 
correlation functions. The same effects, perhaps including some collectivity due to a small 
but finite source motion, may cause the slight decrease towards larger q that is not observed 
in the momentum-integrated correlation functions (cf. Fig. 1). We observe, however, that 
the peak height at g ~ 20 MeV/c relative to the uncorrelated background at g > 40 MeV/c 
is virtually identical in all four cases. This is consistent with what is expected for the ideal 
situation of a purely statistical source, confirming that effects caused by collective radial 
motion or emission times related to the temporal evolution of an expanding source should 
be small for spectator decays. 

The choice of a threshold of Eiab > 20 MeV, on the other hand, raises the question of 
what reaction stage is mainly represented by the obtained correlation data. Light-particle 
spectra from the same experiment, measured with high-resolution telescopes at backward 
angles, give evidence for emission prior to the final breakup stage [^ . This was concluded 
from the comparison of the slope temperatures and multiplicities to the predictions of the 
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Figure 3: Proton-proton corre- 
lation functions, integrated over 
Z bound, for four intervals of the 
absolute value of the proton pair 
momentum. The interval 70 < 
q < 90 MeV/c was used for nor- 
malization. 



statistical multifragmentation model. Distinct components of equilibrium yields and of 
faster particles, termed pre-equilibrium or first stage, have also been identified by other 
groups in data for comparable reactions [O, p3[. In all these cases the main parts of 



the low-energy equilibrium components are below the present threshold. Long-lifetime 
components are thus excluded from the analysis so as to preserve the sensitivity of the 
correlation function to the initial spatial dimensions. On the other hand, pre-equilibrium 
or pre-breakup particles as they are scattered from the forming spectator matter may also 
contribute to an interferometric picture that reflects the extension of the latter. 

Besides the p-p correlations also the p-a and d-a correlations were used to determine 



breakup radii by comparing them to the numerical results of Boal and Shillcock [Q . The 
calculated correlation functions which in their work are given for a discrete set of source 
radii were interpolated, and a Monte Carlo procedure was used for an event-by-event 
simulation of the effects of multiple scattering in the target and of the spatial and energy 
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Figure 4: Breakup density p/po as deduced from the p-p, p-^He, and d-^He correlation 
functions. The rapid decrease of the density with decreasing Zbound reflects the changing 
mass of the spectator source. The error bars represent the statistical uncertainty of the 
radii as given in the table (see text). Note the logarithmic ordinate scale. 



resolution of the detection system p5[. The Coulomb acceleration by the emitting source 
was also here neglected. This seems justified for d-a because the effects should not be 
much larger than in the p-p case. For p-a, the different charge-to- mass ratios of the two 
particles may be expected to cause larger distortions but, at the same time, the finite 
lifetimes of these resonances may reduce them considerably. The resulting fits to the data 
are rather satisfactory (Fig. 1). The corresponding radii are about 1 to 2 fm larger than 
those deduced from the p-p correlations (Table 1). We also observe a slightly stronger 
variation with the pair momentum than in the p-p case. For p-a, e.g., the peak of the 
resonance grows from about 1.3 to 1.6 as the pair momentum is varied from Pgum < 800 
MeV/c to Psum > 1000 MeV/c (result after integration over Zbound)- 

Densities were calculated by dividing the number of spectator constituents, taken from 
the calorimetric analyses of Refs. [^, |2^ and listed in Table 1, by the source volume. The 
densities vary considerably with centrality (Fig. 4) even though the radii are approximately 
constant. This is caused by the varying spectator masses which, in excellent agreement with 
the prediction of the geometric participant-spectator model [^, decrease with increasing 



centrality almost in proportion to Zhound- 

In the p-p case, the mean relative densities decrease from p/po = 0.4 for the near- 
peripheral to below p/po = 0.2 for the most central collisions. These values compare 
well with the densities assumed in the statistical multifragmentation model, including 
their variation with centrality. In the model, the mean density changes as a function of 
the multiplicity if a fixed socalled crack width is used as criterion for the placement of 
fragments inside the breakup volume ||^. For the most peripheral bin, the interpretation 
of the shape of the measured p-p correlations is somewhat uncertain (cf. Fig. 1), and the 
corresponding density has not been plotted. For p-a and d-a, the fits are satisfactory for the 
largest Zf,ound but correspond to rather large radii (Table 1). The densities are therefore low 
which is somewhat unexpected as here the production of highly excited heavy residues is 
the dominant reaction channel |^. Obviously, the assumptions of a homogeneous spherical 
source and of a rapid volume breakup are less well justified in this case which causes 
difficulties for the interpretation of the obtained density values. 

The differences between the results for p-p and for the resonances involving a particles 
are statistically significant but their origin is not clear at present. Besides the systematic 
uncertainties of the employed formalisms, there is also the possibility of differences in the 
scattering cross sections; larger cross sections will cause larger apparent source sizes. An 
interesting connection exists between the large radii derived for the lithium states and the 
limitation of the temperatures deduced from their relative populations [^ . A late emission 
of these resonances could explain both observations. 

In summary, correlation functions constructed from proton-proton coincidences and 
from coincidences between protons, deuterons, or tritons with a particles consistently show 
that the breakup volume does not appreciably change with impact parameter even though 
the spectator mass varies considerably. A quantitative analysis, in the limit of zero lifetime 
of the source, yields results that are consistent with a very moderate radial expansion. The 
deduced breakup densities are rather low as assumed in the statistical model scenarios. The 
variation of the density with impact parameter is caused by the changing spectator mass. 
With regard to the method, it has become clear that the systematic uncertainties of deriving 
densities from correlation functions may be large, in particular for the resonances involving 
a particles for which the formalism is not yet as much advanced. It seems, however, that 
the spectator decay at relativistic energies may represent a particularly favorable case 
because the source charge is moderate, collective motion is nearly nonexistent, and with 
the choice of high energy thresholds long-lifetime components may have been efficiently 
excluded from the analysis. 

The authors would like to thank R. Lednicky for making his code available to them 
and W.A. Friedman for fruitful discussions. M.B., J. P., and CS. acknowledge the financial 
support of the Deutsche Forschungsgemeinschaft under the Contract No. Bel634/1-1, 
Po256/2-l, and Schw510/2-l, respectively. This work was supported by the European 
Community under contract ERBFMGECT950083. 

9 



References 

[a] Present address: Fachbereich Physik, Bergische Universitat, D-42119 Wuppertal, Ger- 
many 

[b] Present address: Dept. of Physics, Yale University, New Haven CT 06512, USA 

[c] Present address: Max-Planck-Institut fiir Kernphysik, D-69117 Heidelberg, Germany 

[d] Present address: Max-Planck-Institut fiir Physik, D-80805 Miinchen, Germany 
[1] D.H.E. Gross, Rep. Prog. Phys. 53, 605 (1990) 

[2] J.P. Bondorf et al, Phys Rep. 257, 133 (1995) 

[3] See, e.g., H. Jaqaman et ai, Phys. Rev. C 27, 2782 (1983) 

[4] For a recent review see L.G. Moretto and G.J. Wozniak, Ann. Rev. Nucl. Part. Science 
43, 379 (1993) 

[5] W. Reisdorf and H.G. Ritter, Ann. Rev. Nucl. Part. Science 47, 663 (1997) 

[6] A. Schiittauf et al, Nucl. Phys. A 607, 457 (1996) 

[7] D.R. Bowman et al, Phys. Rev. Lett. 67, 1527 (1991) 

[8] J. Hubele et al, Phys. Rev. C 46, R1577 (1992) 

[9] K. Hagel et al, Phys. Rev. C 50, 2017 (1994) 

[10] A.S. Botvina et al, Nucl. Phys. A 584, 737 (1995) 

[11] U. Milkau et al, Phys. Rev. C 44, R1242 (1991) 

[12] K. Kwiatkowski et al, Phys. Rev. Lett. 74, 3756 (1995) 

[13] J.A. Hanger et al, Phys. Rev. C 57, 764 (1998) 

[14] M. Begemann-Blaich et a/., Phys. Rev. C 58, 1639 (1998) 

[15] S.P. Avdeyev et ai, Eur. Phys. J. A 3, 75 (1998) 

[16] For a recent review and for references to earlier reviews, see D. Ardouin, Int. J. Mod. 
Phys. E 6, 391 (1997) 

[17] J. Pochodzalla et al, Phys. Rev. C 35, 1695 (1987) 



10 



[18] W.G. Gong et al, Phys. Rev. C 43, 781 (1991) 

[19] G. Wang et al, Phys. Rev. C 57, R2786 (1998) 

[20] S. Pratt and M.B. Tsang, Phys. Rev. C 36, 2390 (1987) 

[21] R. Lednicky et al, preprint SUBATECH-94-19, Nantes 1994; L. Martin et al, NucL 
Phys. A 604, 69 (1996) 

[22] Hongfei Xi et al, Z. Phys. A 359, 397 (1997) 

[23] K. Kwiatkowski et al, Phys. Lett. B 423, 21 (1998) 

[24] D.H. Boal and J.C. Shillcock, Phys. Rev. C 33, 549 (1986) 

[25] C Schwarz et al. Proceedings of the International Workshop XXVII on Gross Prop- 
erties of Nuclei and Nuclear Excitations, Hirschegg, Austria, edited by H. Feldmeier 
et al (GSI, Darmstadt, 1999) p. 168 

[26] J. Pochodzalla et al, Phys. Rev. Lett. 75, 1040 (1995) 

[27] C. GroB, PhD thesis, Universitat Frankfurt, 1998 

[28] J. Gosset et al, Phys. Rev. C 16, 629 (1977) 

[29] V. Serfling et al, Phys. Rev. Lett. 80, 3928 (1998) 



11 



